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Abstract: 

Two recent experimental and theoretical studies (Proc. Natl. Acad. Sci. U. S. A. 2014 , 111, 1282-1287; P/iys. Chem. Chem. Phys. 
2014 , 16, 25942-25949) have addressed the problem of tuning molecular charge and vibrational properties of single molecules 
embedded in nanojunctions. These are molecular characteristics escaping so far to an efficient experimental control in broad 
ranges. Here, we present a general argument demonstrating why, out of various experimental platforms possible, those wherein 
active molecules are asymmetrically coupled to electrodes are to be preferred to those symmetrically coupled for achieving a(n 
almost) complete redox process, and why electrochemical environment has advantages over “dry” setups. This study aims at 
helping to nanofabricate molecular junctions using the most appropriate platforms enabling the broadest possible bias-driven 
control of the redox state and vibrational modes of single molecules linked to electrodes. 
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1 Introduction 

In spite of impressive advances in nanoelectronics,—detailed 
characterization and control of molecular properties under in 
situ conditions continue to remain important challenges for 
fabricating and understanding single-molecule junctions. To 
this aim, more recent studies have emphasized the need to go 
beyond electronic transport.—"— Vibrational properties studied 
via inelastic electron tunneling spectroscopy (lETS)— and sur¬ 
face enhanced Raman spectroscopy (SERS)^i^"— may rep¬ 
resent such valuable piece of information, ideally if they are 
acquired concomitantly with transport measurements. This 
has been demonstrated in a recent joint experimental theoret¬ 
ical study on fullerene-based electromigrated junctions^^ and 
a theoretical study— on viologen-based junctions in electro¬ 
chemical (EC) scanning tunneling microscope (STM) setup, 
which was based on experimental data reported previously.— 
Refs. [3 and [3 indicated that a significant tuning on the 
vibrational frequencies and Raman scattering intensities can 
be obtained via bias-driven changes in the charge of the ac¬ 
tive molecule in a current carrying state; tuning the molecu¬ 
lar charge via applied biases enables to control chemical bond 
strengths and, thence, vibrational properties. 
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Obviously, the broadest control that can be achieved pur¬ 
suing this route corresponds to fully change the average re¬ 
dox state of the molecule by adding an entire electron. To 
be specific, we limit ourselves to n-type (LUMO-mediated) 
conduction (orbital energy offset So > 0), as this is the case 
for the molecular junctions of refs. 13 and [3— Reaching 
this ideal limit was im pos sible in the experiments of ref. [isl 
On the other hand, ref. [ij indicated that an almost perfect re¬ 
duced (n fa 1) state can be reached in electrolytically gated 
junctions.— 

Are the obstacles to obtaining a full reduction in some 
molecular junctions merely of technical nature? This is the 
question that initiated the present study. By comparing perfor¬ 
mances of various nanofabrication platforms and identifying 
the ones (depicted in Fig.[^ and Fig.|^ below) that are most ad¬ 
vantageous for achieving an almost complete reduction, the present 
study aims at helping to design molecular junctions enabling the 
broadest bias-driven control over molecular charge and vibrational 
properties. 


2 Model 

The framework adopted here is provided by the single-level Newns- 
Anderson model, which was discusse d and validated for a va¬ 
riety of molecular junctions,i2i^ including those used in the exist¬ 
ing SERS-transport studies. '^4^ We checked that reorganization ef- 
fect s 14,17,18,22- _^ jjqj qualitatively change the present conclusions. 
This is illustrated by the comparison between Fig.j^and Fig.|4] Fig.[T] 
and Fig. and between panels b and c of Fig. To make the pre- 
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sentation as simple and clear as possible, we will give below only 
formulae wherein reorganization effects are disregarded; details on 
the ensemble averaging needed to include these effects are not given 
here hut can be found elsewhere, 

Eqn l[T), (O, (O and Q, which allow to express the cur¬ 

rent 1 and the LUMO occupancy n, constitute the framework of the 
present discussion (e and h are the elementary charge and Planck’s 
constant, respectively). 
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The LUMO occupancy n of the molecule embedded in a biased 
junction comprises contributions ns^t from the two {s and t) elec¬ 
trode s^l 2 ^ 
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The Newns-Anderson model embodies two asymmetries charac¬ 
terized hy the dimensionless parameters y and d defined in eqn ([3j 
and 0 . The potential profile asymmetry (or voltage division factor) 
y quantifies the hias-induced shift in LUMO energy offset relative to 
the Fermi energy {Ef =0) with respect to the case of unbiased elec¬ 
trodes £0 Lj, denote the couplings of the molecule 

to generic “left”/“suhstrate” (s) and “right”/“tip” (t) electrodes. With¬ 
out loss of generality (because it merely amounts to appropriately la¬ 
bel the electrodes) we assume that ^(uhstrate) is not the weakest cou¬ 
pled electrode to the LUMO (Lf < Lj, 0 < 5 < 1/2) and define the 
source-drain bias V with respect to its potential (V = % — %)■ Us¬ 
ing a potential origin located symmetrically between the electrodes’ 
Fermi energies {—e"Vs.t = EeV /2), —1/2 < y < 1/2. y is positive 
(negative) if the bias shifts the LUMO towards the Fermi level of 
electrode i if). 

Similar to / (Fig.|T]and Fig.0, the LUMO occupancy exhibits 
substantial changes within ranges SV ~ L (which are narrow in typi¬ 
cal non-resonant cases F eg) at the biases (case 1 = 0 ) 

Vr.s = eo/(l/2-y) (8) 

Vr, = ~£o/{l/2 + y) (9) 

where the LUMO and electrodes’ Fermi level become resonant and 
plateaus (saturation) above these biases. The saturation values (sub¬ 
script sat) can be expressed as follows (superscripts ± refer to bias 
polarity) 
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3 Results and discussion 


The impact of the source-drain bias on the molecular charge is de¬ 
picted in Fig.[^and Fig.|4] It is worth noting that the asymmetry in 
the electronic coupling (5 7 ^ 1 / 2 , F^ 7 ^ L^) is physically distinct from 
the asymmetry (y 7 ^ 0 ) in the electric potential drop across a molec¬ 
ular junction.iii2i2^ Their impact on n is completely different. The 
comparison between Fig. and b reveals that y 7 ^ 0 merely breaks 
the bias reversal symmetry (n{V) 7 ^ n{—V)). This is the counterpart 
of the current rectification (/(—V) 7 ^ —I{V), cf. eqn (|T]l and 0). As 
long as the molecule-electrode coupling is symmetric (5=1 / 2 ), the 
plateau values at either bias polarity are equal, = 1 / 2 . 

From the perspective of achieving a full redox process in biased 
molecular junctions, this result for symmetric coupling is disappoint¬ 
ing: whatever high is the source-drain bias which a molecular junc¬ 
tion can withstand (|V| V/.j or |V/_/|), whether the potential profile 
is symmetric (y = 0 ) or asymmetric (y 7 ^ 0 ), the average molecular 
charge cannot exceed half of an electron (nj^, = 1/2 for 8 = 1 / 2 , 
cf. eqn 0). This is a surprising result; intuitively, one may expect a 
full change in the molecular redox state (u = 1 ) at biases correspond¬ 
ing to resonant transport (£ 0 ( 1 ^) ^ )ls{V) or /t, (V)).— Infact, exactly 
on resonance (V = Vj- j or the average molecular charge is even 

smaller: n ~ 1 /4, see eqn l ll2b and l ll2b for 5 = 1/2, F fig). 
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According to eqn dlOb and GD, a molecule in a current-carrying 
state can accommodate an average charge larger that n = 1/2 for 
junctions characterized by asymmetric molecule-electrode couplings 
5 7 ^ 1/2. The stronger the asymmetry, the larger is the charge of 
the average redox state (Fig. 0 . An almost perfect reduced state 
(n ~ 1 ) can be achieved for highly asymmetric molecule-electrode 
couplings (5 <C 1) if the potential of the strongest coupled electrode 
(^) is sufficiently negative (sufficiently large positive biases). In such 
cases (Fj F, and V — V/.i F/e), the Fermi energy of electrode .s 
is sufficiently above the LUMO energy, and the electron transferred 
from this electrode with a high rate (Fj) almost entirely remains on 
the LUMO in cases where the transfer from the LUMO to electrode 
t is inefficient (small F(): n Fs/(Fi-|-Ff) = 1 — 5,$.l (Fig. 0). 
This is impossible in cases of symmetric coupling (F* = FO, wherein 
n Rt Fj/ (Fi -f Ff) = 1 / 2 ; on average, half of the entire electron trans¬ 
ferred from electrode i to the LUMO is further transferred to elec¬ 
trode t in cases wherein the two transfer rates are equal (Fig.0). 

The foregoing analysis demonstrated that an almost complete re¬ 
duction can be obtained (i) for a molecule highly asymmetrically cou¬ 
pled to electrodes and (ii) sufficiently high biases. While the first 
condition (F; <C F?) can be easily satisfied in asymmetrical EC-STM 
(see Sec. 14.2b setups, the second condition seems problematic. 
Fig. m illustrates that the biases needed for an almost perfect reduc¬ 
tion are those where current plateaus occur. 

As experiments on molecular junctions did not routinely report 
currents exhibiting plateaus with increasing bias (we are aware of one 
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exception^), this appears to be an important practical limitation. 

A nearly complete reduction would be possible if the LUMO lied 
below the Fermi level of the strongest coupled electrode (n 1 — 5 
for 5 <C 1 form eqn 13) or below the Fermi levels of both elec¬ 
trodes (Fig. [Hi- This results from eqn 13: «~(1 — 5) + 5 = 1 for 
eo < 0 (LUMO below the Fermi level Ef(= 0) of unbiased elec¬ 
trodes). In the absence of any bias, Eq > 0; the Fermi level lies 
within the HOMO-LUMO gap of a molecule linked to electrodes 
(charge neutrality). However, Eq < 0 becomes possible under elec¬ 
trostatic gating. appropriate gate potential Vc (overpo¬ 

tential in electrochemical language, on which the LUMO energy Eq 
linearly depend s 15,25,26,31 ^ lower the LUMO energy below Ef 
(Fig. [6ji). / — V(j (/ -- Eq) transfer characteristics exhibiting max¬ 
ima, which occur at resonance (Eq ~ 0), can be taken as 

indicating a substantial change in the molecular redox state (from 
n 0 for Eq > 0 to n ~ 1 for Eq < 0). This is illustrated in Fig. |3), 
which emphasizes that it is not the (source-drain) bias V, but rather 
the overpotential that determines the reduction efficiency in electro¬ 
chemical environment. 

4 Specific remarks 

In this section we will address three specific issues related to the re¬ 
sults reported above. 

4.1 Rectification vs. reduction 

To avoid possible misunderstandings, in this subsection we briefly 
discuss some of the present results in the context of those reported 
previously. 

In agreement with previous work (e.g., ref. and Is^l . eqn {D, 
(|3 and ([3 show that asymmetric molecule-electrode couplings alone 
(i.e., Fj 7 ^ F, (or 5 7 ^ 1 / 2 ) and y = 0 ) do not yield rectification: the 
aforementioned equations yield I{—V) = —I{V) for y = 0 irrespec¬ 
tive of the value of 8. 

The fact that a potential prohle asymmetry (y 7 ^ 0) yields current 
rectification has been amply discussed. A very incomplete list in¬ 
cludes ref. fi^l2i3i3^ : notice that the parameters p of ref. I29ll3^ 
and a of ref. [3 correspond to that denoted by y in this paper. In par¬ 
ticular, the present eqn Q coincides with eqn (4b) of ref. 0- This 
is why all conclusions on the current asymmetry (“current rectifica¬ 
tion”) emerging from our eqn {D are not new; they coincides with 
those of the earlier works based on the same formula for the current 
1 = I{V). 

Important new results reported in the present paper are that an ap¬ 
plied bias can yield an asymmetric molecular reduction n{—V) 7 ^ 
n{y) and that current asymmetry I{—V) 7 ^ —f(F) and reduction 
asymmetry n{—V) 7 ^ n{V) are conceptually different. In particular, a 
symmetric current I{—V) = —I{V) (no rectification) does not rule out 
an asymmetric LUMO occupancy n{—V) 7 ^ n{y) (cf. Fig.[3 and cur¬ 
rent rectification /(—F) 7 ^ —/(F) does not rule out symmetric plateau 
values of the reduction degree («(F <C —2 Eo) = n{V ^ leg) = 1/2 
forF, = F). 

Concerning the current rectification a final comment is in order, 
however. Similar to ref. |3 and eqn O, dl and 13 refer to sit¬ 
uations wherein reorganization effects are ignored (i.e., A = 0). As 


recent work demonstrated,.^^ in cases where reorganization effects 
are non-negligible (A 7 ^ 0), the asymmetry F^ 7 ^ F/ yields current 
rectification (/(—F) 7 ^ —/(—F)) even if there is no bias-induced shift 
of the LUMO energy (y = 0); see eqn (13), and (18) to (23) of ref. 

Without intending to be exhaustive (current rectification in not 
our primary aim here), we mention that, in cases wherein F^ 7 ^ F, 
and y = 0 , current rectification may also appear due to charging ef¬ 
fects 

4.2 Symmetric vs. asymmetric contact couplings and ex¬ 
perimental platforms 

The above analysis demonstrated the key role played by the 
molecule-electrode couplings (Fx, x = s,t) in determining the reduc¬ 
tion efficiency. For electrodes with wide flat band structures around 
the Fermi energy, Fj^ can be expressed in terms of the electrode den¬ 
sity of states at the Fermi level p, i and effective transfer integrals fj , 
quantifying the charge transfer between electrodes and the dominant 
molecular orbital (LUMO in the specific case considered here)^ 

rx = Pxfx (14) 

Let us now consider the experiment that succeeded to reveal a 
change in the redox state of a Cgg molecule embedded in a biased 
junction by means of simultaneous SERS-transport measurements in 
an electromigration platform. In view of the high molecular sym¬ 
metry and of the (nearly) symmetric electrodes characterizing elec- 
tromigrated junctions, one can assume Ps ~ Pi and L ^ p. This im¬ 
plies Fj ~ Ff, and we have shown above that in this symmetric case 
(S = 1/2) reduction cannot exceed 50%. Drawing (experimental¬ 
ists’) attention on the limited reduction degree that can be achieved 
even if such electromigrated molecular junctions can be brought into 
a current plateaus regime represents an important aim of the present 
theoretical work.— 

In typical STM or CP-AFM experiments, molecules forming regu¬ 
lar self-assembled monolayers (SAMs) are (covalently) bound to the 
substrate. A difference 7 ^ p, (yielding Fj 7 ^ F, via eqn ([Til l) may 
exist because the substrate surface is typically monocrystalline (e.g., 
Au(lll)) while the tip facet is undefined. Still, it is more probable 
that not the difference Ps 7 ^ p, but rather that of the charge transfer ef¬ 
ficiency (ts 7 ^ tf) determines the contact coupling asymmetry F^ 7 ^ F/. 

Quantifying the asymmetry Fj 7 ^ F^ from transport measurements 
is not straightforward; we have seen above that (nearly) symmetric 
curves /(—F) ~ —/(F), which are measured in numerous junctions 
based on symmetric (and occasionally also asymmetric) molecules, 
do not exclude (highly) asymmetric contact couplings. 

A coupling asymmetry F^ > F, has been concluded in ref. EH after 
a detailed analysis of the transport data in CP-AFM junctions. In CP- 
AFM setups such an asymmetry appears to be plausible because one 
can expect that > t,, given the fact that a stable molecule-tip chemi¬ 
cal bond is hard to imagine; normally, charge transport only occurs by 
applying a loading force at the CP-AFM tip.^^ A similar inequality 
(fj > t,, Fj > FQ can also be expected for STM break junctions; the 
formation of a stable molecule-tip chemical bond is implausible dur¬ 
ing repeated processes of rapidly crashing the tip into and retracting 
it away from the substrate.i-^*'^^ So, although per se the asymmetric 
outlook of the STM and CP-AFM setups does not necessarily imply 
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Fj > Ff, this inequality can be expected in view of the different bond 
strength at the contacts. The fact that the symmetry of the electro¬ 
migration platform and the asymmetry of the STM platform do not 
merely refer to the usual schematic illustrations of these setups (like 
those in Fig.|5]and Fig.| 6 ](, but also have a physical content has been 
recently quantified— 

Out of the various experimental platforms employed, EC-STM se¬ 
tups appears to be the most favorable, enabling almost perfect reduc¬ 
tion. In electrochemical environment (EC-STM), the tip typically 
approaches but does not contact the free end of the molecules. The 
different through-bond vs. through-space charge transfer mechanisms 
at the (EC-)STM substrate and tip reflect themselves in significantly 
different transfer integrals (fj ^ f?) responsible for the high asymme¬ 
try Fj ^ Ff (cf. eqn I14l l): in agreement with this analysis, our recent 
works demonstrated very highly asymmetric molecule-electrode cou¬ 
plings: 8 ~ (ref. |2) and 8 10^^ (ref. fT3i . 

4.3 The Newns-Anderson model 

As noted in Introduction, the present study has been motivated by the 
findings of ref. 13 and [3 which demonstrated that reduction is pos¬ 
sible in molecular junctions in current carrying states and concluded 
that the Newns-Anderson model (single level -l- Lorentzian transmis¬ 
sion) represents an adequate theoretical framework. In favor of the 
the Newns-Anderson model one can still add its ability to excellently 
describe the charge transport by tunneling in a variety of molecular 
junctions— The aforementioned represent a sufficient justifica¬ 
tion for adopting the Newns-Anderson model in the present study. 

Still, for reasons delineated below we believe that our main con¬ 
clusion on the impact of coupling asymmetry (F, 7 ^ Fj) on the reduc¬ 
tion efficiency holds beyond the Newns-Anderson framework. 

Using the expressions of the partial LUMO occupancies Uj and n, 
given in eqn llDl and 0, respectively, eqn 0 can be recast as 


Eqn ( fTTt and Gl} yield the following limiting (plateau) values 
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It is noteworthy that both the current, eqn 03, and the plateau 
LUMO occupancies, eqn 03 and Gl can be entirely expressed in 
terms of the rates Fj^. Parenthetically, these equations also holds in 
the presence of reorganization (A 7 ^ 0). In view of the appealing sim¬ 
plicity of the above expressions and their clear physical content, we 
believe that they apply beyond the Newns-Anderson framework, as 
also suggested by the cartoons presented in Fig.|^an dm 

Obviously, the aforementioned should not be taking as attempting 
to discourage alternative approaches of redox processes in molecular 
junctions based on other theoretical models utilized in the literaturei^ 
or by further refinements of the Newns-Anderson model itself. Con¬ 
cerning the latter possibility, an extension that appears to us as par¬ 
ticularly important in studying charge transport through redox units 
is to consider population dependent contact couplings F^, = F^ ,()i) 
and LUMO energy eg = Eo{n). 

From numerous studies on redox electrochemical systems,.?^ it is 
known that changes in the reduction (or oxidation) degree — with 
accompanying redistributions of the electronic charge over the whole 
molecule — have an overall, although selective impact on the chem¬ 
ical bond strengths between the various molecular constituents and 
molecular orbital energies. For molecular junctions, the influence of 
the reduction (or oxidation) degree on the bond strengths at the con¬ 
tacts is of particular interest; so, the transfer integrals , (and thence 
Fj 7 , cf., eqn GUl) are expected to be affected. 


I = 


tts.t — 



(15) 
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In eqn dlSt (where the factor 2 is the spin contribution), the first term 
refers to an electron arriving at the LUMO as a result of the coupling 
of the molecule to the left electrode. Eqn 0 and 0 yield 


F, 1 eo(F)-eF/2 

Us = -arccot-== 1 . 

2F ;r F 


(17) 


This electron is transferred to the right electrode within a character¬ 
istic time T/. Via eqn ( 116b . T, corresponds to the rate Ff determined 
by the coupling to the right electrode. Likewise, the second term of 
eqn (Ell describes the electron flow in opposite direction: an electron 
arriving at the LUMO as a result of the coupling of the molecule to 
the right electrode. Eqn 0 and 0 yield 


n, 


F, 1 ^eo{V)+eV/2 ^ 

-arccot- 0 = 1 1 

IFn F 


(18) 


This electron is transferred to the left electrode within a characteristic 
time Is- Via eqn ( 116b . Tj corresponds to the rate Fj determined by the 
coupling to the left electrode. 


5 Conclusion 

We believe that the new theoretical results reported in the present pa¬ 
per offer relevant information for a better theoretical understanding 
of the microscopic processes occurring in molecular junctions away 
from equilibrium and useful hints for future experimental investiga¬ 
tions. 

We have identified experimental platforms allowing almost com¬ 
plete bias-controlled redox processes in molecular junctions. Infor¬ 
mation on the bias (V and Vg ~ Go) dependent LUMO (in the specific 
case examined) occupancy n, the quantity on which we have focused 
our attention, can be obtained from bias dependencies of vibrational 
frequencies cOy {V, eg) = [1 — n (F, Eg)] cOy + n (F, Eg) COy and Raman 
scattering intensities Ay (F, Eg) = [1 — n (F, Eg)] A" -f n (F, Eg) Ay ex¬ 
tracted from simultaneous SERS-transport measurements, which 
take values interpolating between the relevant (neutral n and an¬ 
ionic a) charge species.— We have demonstrated that the experimen¬ 
tal setup asymmetry plays an essential role in achieving an almost 
perfect reduced state via bias tuning; efficient reduction is possible 
within some molecular electronics platforms, but definitely impos¬ 
sible under other platforms. We have shown that reduction can at 
most reach 50% in two-terminal setups wherein a molecule is sym¬ 
metrically coupled to electrodes (Fig.lSj)). On the contrary, an almost 
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complete redox process can be obtained in cases of highly asym¬ 
metric molecule-electrode couplings (Fig.lSj)). On this basis, an im¬ 
proved reduction can be expected in fullerene-based junctions if an 
STM platform (highly asymmetric coupling) is adopted instead of 
the electromigration platform (symmetric couplings) utilized in ex¬ 
periments.— 

Further, we found that an almost full reduction in two-terminal 
setups is accompanied hy current plateaus. Because such plateau ef¬ 
fects have been observed in molecular junctions based on the zwit- 
terionic molecule hexadecylquinolinium tricyanoquinodimethanide 
(CieFljaQ— 3CNQ)iS (seemingly the only known example), concur¬ 
rent SERS-transport measurements on this system could be of 
interest to investigate bias-driven changes in molecular vibrational 
properties. 

The fact that, unlike existing orbital gating measurements using 
“dry” platforms,^ experiments resorting to electrolyte gating 
succeeded to reveal such a maximum renders the electrochemical 
three-terminal (EC-STM) platform a promising route in achieving an 
efficient reduction in molecular electronic devices,^! also because it 
does not require high biases hardly accessible experimentally and the 
coupling asymmetry is not critical {cf. Fig.|6]5). 
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Fig. 1 Results for eo = 0.5 eV, F = 20 meV, S = 0.1, and 7=0 
showing that the biases V where steps and plateaus in the current I 
and LUMO occupancy n occur are controlled by the values Fj ^ and 
Vt,r (see eqn (5) and (6) of the main text) where the LUMO becomes 
resonant with the Fermi levels of electrodes i(ubstrate) and t(ip), 
respectively. Notice that reduction efficiency (quantified by the 
LUMO occupancy n) is strongly dependent on the bias polarity in 
case of an asymmetric molecule electrode coupling (F^ F,), unlike 

the current plateau value, which is unchanged upon bias polarity 
reversal, and n, represent the separate electrodes’ contributions to 
the LUMO occupancy n expressed by the first and second terms in 
the RFIS of eqn (2) from the main text. For better readability, a 
vertical dashed line marks the zero bias reference. 
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Fig. 2 Results similar to those of Fig.[T]including in addition 
reorganization effects quantified by the reorganization energy X 
whose value is given in the legend. For better readability, a vertical 
dashed line marks the zero bias reference. 




Fig. 3 The stronger the asymmetry of molecule-electrode couplings 
(smaller 5’s), the more efficient is the LUMO reduction. Notice that 
a nonvanishing value of y (panels b and c) breaks the bias reversal 
symmetry (n{—V) ^ n{V)) but does not change the plateau values of 
n with respect to the case 7=0 (panel a). The values of Eq and F are 
comparable to those used in existing concurrent SERS-transport 
studies. For better readability, a vertical dashed line marks the 
zero bias reference. 
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Fig. 5 An almost complete reduction can be achieved if the LUMO 
is much stronger coupled to one electrode (L, S> L/) having a 
sufficiently negative potential (panel a), which is impossible for 
symmetric coupling to electrodes (panel b). 


Fig. 4 Results similar to those of Fig.[3including in addition 
reorganization effects quantified by the reorganization energy X 
whose value is given in the legend. For better readability, a vertical 
dashed line marks the zero bias reference. 
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(a) almost full reduction 
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Fig. 6 Results demonstrating that an efficient reduction can be 
achieved by lowering the LUMO energy below the electrodes’ 
Fermi level via LUMO gating (panel a). An almost complete 
reduction (values of n close to unity) can be obtained for reasonable 
gate-driven LUMO shifts; in this case, the coupling asymmetry 5 
plays a secondary role: n 95% (almost 5-independent) at 
Sq ~ —0.1 eV (without reorganization effects: panel b; with 
reorganization effects: panel c). These results emphasize that it is 
not the (source-drain) bias V, but rather the overpotential (on which 
Eq linearly depends) that determines the LUMO occupancy n in 
electrochemical environment. 
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